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Abstract 

In this paper, we propose two new models in f(T) gravity to realize the crossing of the 
phantom divide line for the effective equation of state, and we then study the observational 
constraints on the model parameters. The best fit results suggest that the observations favor a 
crossing of the phantom divide line. 
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I. INTRODUCTION 

Various observations jl-4| have confirmed the fact that our Universe is undergoing an 
accelerating expansion and it entered this accelerating phase only at the near past. The 
proposals that have been put forth to explain this observed phenomenon can basically 
be classified into two categories. One is to assume the existence of an exotic energy 
with negative pressure, named dark energy. The simplest candidate of dark energy is the 
cosmological constant with the equation of state w = — 1 |5j. It, however, suffers from two 
serious theoretical problems, i.e., the cosmological constant problem and the coincidence 
problem. Thus, some scalar field models, such as quintessence j^] and phantom [7], are 
proposed. For single scalar field models, it has been shown that the equation of state 
cannot cross the phantom divide line (w = —1). So, models with a combination of 
phantom and quintessence 8 , [t| , and scalar field models with scalar-dependent coupling 



1 Of ] as well as fluid models [ri| have also been constructed to 



in front of kinetic term 
realize the crossing of the phantom divide line, which still seems to be allowed by recent 



observations 



12 
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Another alternative to account for the current accelerating cosmic expansion is to 

n 

modify Einstein's general relativity theory. One such modification is the f{R) theory [14] 
(see 15] for recent reviews), where the Ricci scalar R in the Einstein-Hilbert action is 
generalized to an arbitrary function / of R. For this theory, it has been found that 
the effective equation of state can cross the phantom divide line from phantom phase to 
non-phantom one 16|, ll7 ] . 

Recently, a new modified gravity which can also explain the accelerating cosmic ex- 
pansion 18], named f(T) theory, has spurred an increasing deal of interest. The f(T) 
theory is obtained by extending the action of teleparallel gravity [19| in analogy to the 
f(R) theory, where T is the torsion scalar. An important advantage of the f(T) the- 
ory is that its field equations are second order as opposed to the fourth order equations 

rn 

of f(R) gravity. More recently, Linder [20|] proposed some concrete f(T) models (see 



also Ref. 



21|1. We placed observational constraints on the parameters of some of these 



models [20, |22j, [23(, in particular, and analyzed the dynamical properties of the f(T) 
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theory [24j ]. in general. A reconstruction of the f(T) theory from the background ex- 
pansion history and the /(T) theory driven by scalar fields were studied in 25], and the 
cosmological perturbations and growth factor of matter perturbations in the f(T) theory 



were investigated in Refs. 
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examined in Re 
Refs. 
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,122| 
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s. 
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271 ] . In addition, the issue of local Lorentz invariance was 
30]. It should be noted, however, that the analysis performed in 



3jJ indicate that models proposed so far in the f(T) theory [20j,]2l|] behave 
quintessence-like or phantom-like, and thus cannot realize the crossing of the phantom 
divide line for the effective equation state, although the observational data 
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seems 



to indicate this crossing is still a possibility not ruled out. So, in this paper, we propose 
two new f(T) models which can realize the crossing of —1 line, and then discuss the 
constraints on model parameters from recent observations. A remarkable feature of our 
models is that they realize the crossing of the phantom divide line from a non-phantom 
phase to a phantom phase in contrast to the viable f(R) models where the phantom 



divide line is crossed the other way around 



171 ] . It is interesting to note that a crossing of 



the phantom divide from the non-phantom phase to the phantom one is consistent with 
the recent cosmological observational data 12j. Finally, let us note that, recently, a new 
model with the crossing of phantom divide line is also proposed in 



II. THE f(T) THEORY 

The f(T) theory is obtained by extending the action of teleparallel gravity to T + 
f(T). The teleparallel theory of gravity is built on teleparallel geometry, which uses the 
Weitzenbock connection rather than the Levi-Civita connection. So, the spacetime has 
only torsion and is thus curvature-free. 

Assuming that the universe is described by a flat homogeneous and isotropic 
Friedmann- Rob ert son- Walker metric 

= diag(l, -a 2 (t), -a 2 {t), -a 2 (t)) , (1) 

where a is the scale factor, it has been found in Refs. 
the teleparallel gravity can be expressed as 

T = -QH 2 , (2) 
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341 ] that the torsion scalar in 



with H = da 1 being the Hubble parameter. In addition, the modified Friedman equations 
have the following form 

2 8nG f T 
#=^p-g+3/T, (3) 

2 16rrGP -T + f-2Tf T 
{H)= -4Tf TT -2-2f T ' (4) 

where a prime denotes a derivative with respect to In a, the subscript T represents a 
derivative with respect to T, p is the energy density and P is the pressure. Here we 
assume that there are both matter and radiation components in the Universe, thus 

P = Pm + Pr, P = T^Pr • (5) 

If we rewrite the modified Friedmann equation (Eq. (j3j)) in the standard form as that 
in general relativity, we can define an effective dark energy, whose energy density can be 
expressed as, 

Here > 1 is required in order to have a positive value for p e jf- This usually gives a 
constraint on physically meaningful models. Using energy conservation equation, p e ff + 
3H(1 + w e ff)p e ff = 0, one can yield the effective equation of state w e ff 

f/T -f T + 2Tf TT + IjgHfr + 2Tf TT ) 
Weff (l + f T + 2Tf TT )(f/T-2f T ) ■ 1 j 

The same expression can also be obtained using Eq. (j3J) to define an effective pressure 
p e ff and then deriving w e ff. 



III. TWO NEW f(T) MODELS 

In this section, we propose two new f(T) models, labeled as Model A and Model B, 
which can realize the crossing of the phantom divide line for the effective equation of 
state. 
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• Model A 

/(T) = a (-T)«tanh^, (8) 

where a and n are two model parameters. The requirement of 2 -y x > 1, which ensures 
p e ff > 0, gives rise to 2n — Ax ■ csch[2x] > 1 with x = Since 4x • csc/i[2a;] < 2, n must 
be greater than |. Substituting Eq. (|H]) into the modified Friedmann equation, we have 

a ~ ~(6i/ 2 )"- 1 [2sech(l) 2 + (1 - 2n) tanh(l)] ' 1 J 

Here Q m0 and f2 r0 are the present dimensionless density parameters of matter and radia- 
tion, respectively. 

In Fig. (JT]), we show the evolutionary curves of the effective equation state with different 
values of n (right panel) and the cosmic evolution with n = 1.65 (left panel). From the 
right panel, one can see that the effective equation of state firstly crosses the phantom 
divide line from > — 1 (non-phantom phase) to < — 1 (phantom phase), and then evolves 
to > —1. So, it crosses the —1 line twice. In order to illustrate why this phenomenon 
occurs, we plot a figure (FigfSJ) to show the regions w e ff < — 1 in n — E plan with 
Qmo — 0.26, where E = H/Hq. From this figure, one can see that n must be smaller 
than a critical value, i.e. n < 1.686 when Q mQ = 0.26, to render w e ff cross —1, and, 
once w e ff cross the —1 line, it must cross it twice. This makes the f(T) models distinct 
from the viable f(R) models where only a crossing from phantom phase to non-phantom 



one is allowed 



17| . Finally, w e // approaches to —1, which means that the final state 



of our Universe is an exponential expansion phase. This result is consistent with what 



obtained in Ref. 24| where it has been found through the dynamical analysis that the 
Universe in the f(T) theory finally enters a de Sitter expansion phase. Furthermore, 
the right panel reveals that the Universe, in this model, has a long enough period of 
radiation domination to give the correct primordial nucleosynthesis and radiation-matter 
equality, and a matter dominated phase. In other words, the usual early universe behavior 
can be successfully obtained to agree with the primordial nucleosynthesis and the cosmic 
microwave background constraints. 
• Model B 

f(T) = a(-T) n (l - e pTo/T ) (10) 



z 



ln(l+z) 



FIG. 1: The evolutionary curves of the effective equation of state with different values of n 
and il m o = 0.26 (left panel), and the cosmic evolution with n = 1.65, Vl m Q = 0.26 and il r o = 
0.26/1200 (right panel) for Model A. In the right panel, the dot-dashed, dashed, and solid lines 
represent the evolutionary curves of the dimensionless density parameters for the effective dark 
energy, radiation and matter, respectively. 
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FIG. 2: The regions of phantom and quintessence in n — E plane with Q m o = 0.26, where 
E = H/Hq. The red dashed line is the minimum value to which the universe can reach. 

with three model parameters a, n and p. From ^j 1 - > 1 given by the requirement of 
p e ff > 0, we obtain that 2n — jjrj^ > 1 with x = pT /T. This leads to n > 0.5 since 
— 2 *+ eX < 0. So, we now restrict our discussion to the case of n > \ for Model B. Using 




6 



the modified Friedmann equation, we have 

(6#*)i-«(i _ n m0 - a r0 ) 



a 



-1 + 2n + eP(l - 2n + 2p) 



The exponential model given by Linder 20J is a special case, the n = 1 case to be exact, of 
the present model. When p — 0, our model reduces to the power low model f(T) ~ (— T) n , 



18 



20 



22 



24(. Let us note that when n = 1 



which has been studied in detail in Refs. 
or p = 0, the crossing of phantom divide line is impossible as we have already pointed 



out |24J |. This is also confirmed by the n = 1 case in Fig. ([3]) in the present paper. Fig. Q 
shows the evolutionary curves of w e ff (left and middle panels) and the comic evolution 
(right panel) for model B. We find that, for the crossing of the —1 line to occur, it is 
required that p and n — 1 should have the same sign. When p > and n > 1, w e ff evolves 
from > —1 to < —1, while, when p < and | < n < 1, the crossing direction is just the 
opposite. In addition, we also find that, the model behaves like quintessence when p > 
and | < n < 1, and like phantom when p < and n > 1. The right panel in Fig. (|3j) 
gives the comic evolution with n — 1.1 and p = 0.1, from which one can see that the usual 
early universe behavior can also be obtained just as Model A. 
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FIG. 3: The evolutionary curves of the effective equation of state with different values of p and 
n and f2 m o = 0.26 (left and middle panels), and the cosmic evolution with n = 1.1, p = 0.1 
fi m o = 0.26 and f2 r o = 0.26/1200(right panel) for Model B. In the right panel, the dot-dashed, 
dashed, and solid lines represent the evolutionary curves of the dimensionless density parameters 
for the effective dark energy, radiation and matter, respectively. 
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IV. OBSERVATIONAL CONSTRAINTS 



Now, we discuss the constraints on model parameters of Model A and Model B from 
recent observational data, including the Type la supernovae (Sne la), the baryonic acoustic 
oscillation (BAO) distance ratio and the cosmic microwave background (CMB) radiation. 
The Sne la data used in our analysis is the Union2 compilation released by the Supernova 



Cosmology Project collaboration recently 35|, which consists of 557 data points and is 



the largest published sample today Using the usual method, we constrain the theoretical 
model from the Sne la by minimizing the y 2 value 

-2 _ \^ \Pobs(Zj) - fith{Zj)} 2 , v 

XSne - ' y LZ > 

i=l «>* 

where o 1 i are the errors due to the flux uncertainties, intrinsic dispersion of Sne la absolute 
magnitude and peculiar velocity dispersion. fi obs is the observed distance moduli and \i t h 
is the corresponding theoretical one, which is defined as 

H th = 5\og w D L - no . (13) 

Here fiQ = 5 log 10 /i+42.38 with h = HQ/100km/s/Mpc, and Dl is the luminosity distance, 

D * s < 1+ *£m- (14) 

with E(z) = H(z)/Hq. In order to marginalize the nuisance parameter /i (or h), following 



the approach given in Ref. |36|], we expand x 2 Sne to xIne(A*o) = — 2Bfi + C with 
A = E V< i; B = Y.[^obs(zi) - 51og 10J D L ]/<7^ and C = J2[Vobs(zi) - 5 log 10 D L f/a 2 u ^ 
and find that x| ne has a minimum value at [1q = B/A, which is given by 

Xlne = C ~J[- ( 15 ) 

Thus, we can minimize Xsne instead of Xs ne to obtain constraints from Sne la. 

For the BAO data, the BAO distance ratio at z = 0.20 and z = 0.35 from the joint 



analysis of the 2dF Galaxy Redsihft Survey and SDSS data 33] is used. This distance 
ratio 

D v (z = 0.35) 



D v (z = 0.20) 



1.736 ±0.065 (16) 



is a relatively model independent quantity with D v (z) defined as 



D v {z b ao) 



zbao ( f ZBAO dz ^ 1/3 



(17) 



H(zbao)\Jo h ( z 
So, the constraint from BAO can be obtained by performing the following x 2 statistics 

2 [Dy(z = 0.35)/D v (z = 0.20) - 1.736] 2 

Xbao ~ • (18) 



Finally, we add the CMB data in our analysis. Since the CMB shift parameter R |38l.l39] 
contains the main information of the observations from the CMB, it is used to constrain 
the theoretical models by minimizing 



2 [R - R obs } 2 

XCMB - ~2 > l iy J 



0, 



where R Q b s = 1.725±0.018 [40J, which is given by the WMAP7 data, and its corresponding 
theoretical value is defined as 

fZCMB A-l 

*=<&l 4k' (20) 

with zcmb = 1091.3. 

The constraints from a combination of Sne la, BAO and CMB can be obtained by 
calculating x| ne + Xbao + Xcmb- We find that, for Model A, the best fit values occur at 
Q m0 = 0.282 and n — 1.65 with Xmn = 543.948. The contour diagrams at the 68.3% and 
95.4% confidence levels are given in Fig. (jlj). From this figure and Fig. <§Q), we conclude 
that the observation favors a crossing of phantom divide line. 

For Model B, the best fit values of model parameters are Q m o = 0.267, p = 0.02 and 
n = 1.08 with x\n n = 544.213. It is easy to see that the best fit value favors a crossing 
of the phantom divide line from > — 1 (non-phantom phase) to < — 1 (phantom phase). 
This is consistent with the recent observational data [12J but is opposite to what was 
found in viable f(R) models [171 ] . Fig. ([S]) gives the constraints in the n—p plane with 
Q m0 = 0.267 at the 68.3% and 95.4% confidence levels, and in this figure n > \ given by 
the requirement of p e ff > has been taken into consideration. From Figs. (J3J , one can 
see that all possible behaviors of w e ff shown in the left and middle panels of Fig. (j3]) are 
allowed by observations. 

9 
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FIG. 4: The constraint on £l m o and n at the 68.3% and 95.4% confidence levels for Model A 
from Sne Ia+BAO+CMB. 

Now we consider the constraints on the ACDM model. The best fit result is Q m o = 
0.270 with xiiin = 544.403. This x\ti n * s slightly larger than that obtained in the above 
two f(T) models. With the x\ lin jdof (dof: degree of freedom) criterion, the ACDM is 
slightly favored by observations. 
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FIG. 5: The constraint on p and n for Model B with ft m0 = 0.267 at the 68.3% and 95.4% 
confidence levels from Sne Ia+BAO+CMB. n > ^ given by the requirement of p e ff > is 
considered. 
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V. CONCLUSION 



The f(T) theory is a new modified gravity, obtained by extending the teleparallel 
gravity, to account for the current accelerating cosmic expansion without the need of 
dark energy. In this paper, we have proposed two new f(T) models in which the crossing 
of the phantom divide line is possible. A remarkable feature of the our models is that they 
realize the crossing of the phantom divide line from a non-phantom phase to a phantom 
phase in contrast to the viable f(R) models where the phantom divide line is crossed 
the other way around |l7j. It is interesting to note that a crossing of the phantom divide 
from the non-phantom phase to the phantom one is consistent with the recent cosmological 



observational data 



121 ] . By studying the evolutionary curves of w e ff, we find that w 



eff 



can cross the —1 line in both models and it is crossed twice in Model A. Furthermore, we 
also find that both models can produce the usual early universe behaviors in the sense that 
they both allow a long enough period of radiation domination and a matter dominated 
phase to agree with the primordial nucleosynthesis and the cosmic microwave background 
constraints. We have also discussed the constraints on model parameters from recent 
observations including Sne la, BAO and CMB. Our results show that observations favor 
a crossing of the —1 line for Model A, whereas, for Model B, all possible evolutions for 
w e ff given in Fig. (j3J) are allowed, although the best fit result favors the crossing. With 
the XMin/dof (dof: degree of freedom) criterion, we find that the ACDM is still favored 
slightly by observations. 
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